Controlling non-Abelian statistics of Majorana fermions in semiconductor nano wires 
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Under appropriate external conditions a semiconductor nanowire in proximity to an s-wave superconduc- 
tor can be in a topological superconducting (TS) phase. This phase supports localized zero-energy Majorana 
fermions at the ends of the wire. However, the non-Abelian exchange statistics of Majorana fermions is difficult 
verify because of the one-dimensional topology of such wires. In this paper we propose a scheme to transport 
Majorana fermions between the ends of different wires using tunneling, which is shown to be controllable by 
gate voltages. Such tunneling-generated hops of Majorana fermions can be used to exchange the Majorana 
fermions. The exchange process thus obtained is described by a non-Abelian braid operator that is uniquely 
determined by the well-controlled microscopic tunneling parameters. 
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I. INTRODUCTION 

Majorana fermions (MFs) have been the subject of in- 
tense recent study in part due to their potential application in 
topological quantum computation (TQC). 17 Unlike ordinary 
fermionic or bosonic operators for which the particle creation 
operators are the hermitian conjugate of the annihilation op- 
erators, MF operators, 7, are self-adjoint (7^ = 7). In this 
sense, MFs are their own anti-particle and the realization of 
such excitations would be the first example of such particles 
which had been proposed more than seventy-five years ago. 8 
MFs are of interest for TQC because despite having no inter- 
nal degrees of freedom individually, a pair of MFs, say 71 and 
72, have two distinct possible states (fusion channels). These 
states, which may be thought of as the two possible occupa- 
tion states of the complex fermionic operator c 1 = 71 + 2 n2 , are 
energetically degenerate to a degree exponential in the sepa- 
ration of the Majorana fermions, and correspond to the eigen- 
states of the combined operator 17172, with eigenvalues ±1. 
Such a topological protected degeneracy has yet to be seen in 
nature, and the observation of such would be a major break- 
through in physics. The central idea of TQC is to use the 
topologically degenerate states of a pair of MFs as a 2-level 
topological qubit which would in principle be protected from 
decoherence. The manipulation of the information contained 
in the topological qubits requires the use of topological braid 
operations which consists of moving the MFs around one an- 
other. 

In the past few years, topological superconductors have 
become promising candidates for realizing MFsi2~— Re- 
cently, it has been proposed that a semiconductor thin film 
with Rashba-type spin-orbit (SO) coupling together with 
proximity-induced superconductivity and Zeeman splitting 
would be a suitable platform for realizing a Majorana- 
fermion-carrying topological superconducting (TS) stated— 
The s-wave superconducting pairing potential can be induced 
in the semiconductor system by placing it in proximity to a 
conventional superconductor such as aluminum. The Zee- 
man splitting can similarly be induced, in principle, by prox- 
imity to a magnetic-insulator. 20 The one-dimensional ver- 



sion of this system, i.e. a semiconducting nanowire with 
proximity-induced s-wave superconductivity, has also been 
shown to host MF as zero-energy modes at the ends of the 
wire under appropriate conditions. 21,22 The one-dimensional 
nanowire geometry has the specific advantage that the Zeeman 
splitting Vz in the nanowire is not required to be proximity- 
induced from a magnetic insulator, but instead can be intro- 
duced by a magnetic field parallel to the nanowire. 22 Such a 
parallel magnetic field would not introduce unwanted orbital 
effects such as vortices if a thin-film superconductor is used 
to generate the superconducting proximity-effect. The pro- 
posed semiconducting structures exist in a TS phase and sup- 
ports MFs at its ends when the s-wave superconducting pair 
potential A, Zeeman splitting Vz, and the chemical poten- 
tial fj, satisfy the condition V| > A 2 + ^ 2 . 17 i 20 i 22 Thus the 
chemical potential /i, which can be controlled by an external 
gate potential, can be used to tune a nanowire from the TS 
phase to a non-topological (NTS) phase. In fact, the s-wave 
proximity effect on a InAs quantum wire, which also has a 
sizable SO coupling, may have already been demonstrated in 
experiments. 23 These semiconductor based proposals for real- 
izing a TQC platform can take advantage of the considerably 
advanced semiconductor fabrication technology. Therefore, it 
seems that a Majorana-carrying TS state in a semiconductor 
quantum wire may be within experimental reach. 

Until recently, motivated by experiments, most discussions 
of observing non-Abelian statistics using MFs has been re- 
stricted to 2D systems. In 2D non-Abelian systems, the quan- 
tum information associated with MFs can be manipulated in 
a topologically protected manner by exchanging the Majo- 
rana bound states (e.g., by adiabatically moving vortices in 
p+ip superconductors) i^i 2 ^ The protection of the topological 
degeneracy associated with MFs requires the MFs to remain 
spatially separated at all steps of the exchange. Therefore, at 
first glance, it appears that it is impossible to exchange the 
MFs at the ends of a ID wire, since any such attempt would 
necessarily lead to the spatial overlap of MFs at some stage 
of the overlap process. A solution to this problem has been 
provided by by Alicea et al2&, who have shown that connect- 
ing up a system of nanowires into a network allows one to 
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exchange MFs. In their proposal, the two MFs, 71 and 72, at 
the ends of a wire can be exchanged by introducing an addi- 
tional nanowire B. The additional nanowire, B, allows one 
to temporarily move one of the MFs, say 71, away from the 
original nanowire, A, so that the other MF 72 can be moved 
across the wire A without colliding with 72. The MF 71 can 
then be returned from the wire B back to it's original wire A. 
While this scheme solves the basic problem of non-Abelian 
statistics in ID, it requires the transport of a MF across a tri- 
junction between two topological nanowires A and B, which 
is potentially a more complex topological object than the sim- 
ple topological nanowire. The continuous transport of MFs 
through such a tri-junction is potentially dependent on details 
of the junction that may be difficult to control. 27 Moreover, 
from a theoretical point of view, the explicit determination of 
the non-Abelian statistics in this geometry in terms of the mi- 
croscopies of the junction is somewhat complicated. 

In this paper we propose an alternative scheme to transport 
MFs at the ends of one dimensional semiconductor nanowires 
where the ends of the nanowires remain fixed, but the tunnel- 
ing amplitudes between the end MFs are varied. Bringing the 
end MFs closer together allows one to create a non-vanishing 
Hamiltonian of the MFs which can generate effecting MF 
hopping from one end site to another. Using this picture of 
a dynamically changing tunneling Hamiltonian, we will be 
able to derive a simple explicit expression for the non-Abelian 
statistics transformation of the MFs in terms of tunneling ma- 
trix elements. 



II. OUTLINE AND SUMMARY OF RESULTS 

As mentioned in the previous paragraph, MFs are strictly 
zero-energy modes with an associated topological degeneracy 
only in the limit when they are separated by a distance that 
is large compared to the decay length of the MFs. The two 
states of the MFs 71 and 72 can be described in terms of the 
2 possible occupation states of the Dirac fermion c) = 71 + 
«72- These 2 states correspond to the eigenvalues and 1 of 
the number operator n = c)c = 1+I ^ 172 . In general, the 
Hamiltonian for a pair of MFs with a non-negligible splitting 
produces a splitting between the 2 energy states and can be 
written as 

Htunneling = i(l2(x)jl"/2 (1) 

where (i2(x) is the tunneling matrix element for the MFs 
which depends on the separation x between the MFs 71 and 
72^ The energy splitting between the n = and n = 1 states 
is given by | Ci ,2 (a^) I ■ Therefore, the topological degeneracy 
of MFs emerges only in the limit x £ when the MF overlap 
matrix-element £12 (x) vanishes because of the localization of 
the MF wave-functions. Here £ is the localization length of 
the MFs. The tunneling of MFs at the ends of different wires, 
whose ends are placed close together, is entirely analogous to 
the tunneling of electrons between two quantum dots which 
can be controlled by raising and lowering the barrier between 
the dots. Similarly, tunneling amplitudes between MFs on dif- 
ferent semiconductor nanowires can be controlled simply by 



adding a gate controllable tunnel barrier between the MFs. 
Gate voltages can also induce tunneling between MFs at the 
ends of the same TS segment by tuning the nanowire close to 
a TS-NTS phase transition. 20 Bringing the nanowire close to 
the TS-NTS transition, decreases the gap of the system which 
in turn increases the localization length £ of MFs in the wire 
and allows the tunneling between the initially localized MFs 
at the ends of TS segments. The quantitative details of how the 
tunneling is controlled in topological nanowires is discussed 
in the appendix. 

Tunneling of ordinary fermions such as electrons can be 
used to move electrons from one quantum dot to another in 
a system of quantum dots. In this paper, we will show that 
the same principle applies to MFs, and repeated use of the 
tunneling Hamiltonian in Eq.[T]can be used to exchange MFs 

71 and 72 in a system of TS nanowires that hosts such MFs 
at its ends. The unitary time-evolution operator U associ- 
ated with the exchange maps 71 — > UjiW = A72 and 

72 — > U^zU' = — A71 where A can be directly computed 
from the tunneling matrix elements Qj involved in moving 
the MF 71 to the starting position of the MF 72. In the low- 
energy subspace of MFs, the time-evolution operation U, that 
describes the exchange process has the usual form of a braid 
matrix 2 ^ 

U = e^ Xl112 . (2) 

While there are only two possible answers A = ±1 for the 
braid-matrix, it is critical to be able to determine the factor 
A for a given braid since this is what distinguishes 'clock- 
wise' from 'counter-clock- wise' exchanges. For the specific 
geometry discussed in the appendix with wires placed in a 
superconducting film together with an in-plane magnetic field 
at 45 degrees to the wires, the sign of the braid matrix, A, 
is determined by the sign of the Rashba spin-orbit coupling 
constant a. 

The braiding scheme we will discuss is potentially re- 
lated to measurement-only schemes for braiding of topolog- 
ical quasiparticlesi 2 ^ 2 ^ However, it is not clear how the anyon 
model postulates assumed in the measurement-only theory ap- 
ply to the superconducting nanowire systems described by 
mean-field BCS theory. For example, the identification of 
the tunneling matrix element between MFs in Eq. Q]with the 
topological charge measurement in Ref. 29 becomes subtle 
in cases where the sign of the tunneling £12 (aO oscillates in 
sign with the separation x. On the other hand the approach 
in this paper is based only on the MF tunneling Hamilto- 
nian Eq. Q] which can be derived microscopically from BCS 
Hamiltonians. 31 The tunneling matrix elements Cij(x) them- 
selves depend on the details of the nanowire system such as 
the spin-orbit coupling, the orientation of the wire and the 
Zeeman splitting. Therefore, we first consider exchange of 
MFs around a specific triangular loop geometry in terms of the 
tunnel matrix elements between the various MFs and then in 
the appendix, we show how the microscopic tunneling param- 
eters may be calculated in one specific geometry. As a result 
of our calculation, we find that for general values of the tun- 
neling, the parameter A in the braid matrix U has the simple 
form A = sgn((i2)x where \ is the junction chirality of the 
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FIG. 1 : (Color online)(a) Combination of SC, Zeeman and gate po- 
tentials leads to nanowire segments in TS and NTS phases. The 
gates 1 and 3 are adjusted such that the nanowire is in TS phase, 
while 2 is adjusted so that the wire is in NTS phase corresponding to 
the schematic in panel (b). (b) Nanowire segments in TS phase are 
shown as blue (solid) lines and NTS segments are shown as red (dot- 
ted) lines. Orange (light) and blue (dark) circles indicate unpaired 
and paired MFs at TS-NTS interfaces, respectively. MFs are paired 
by tunneling across the TS or NTS segments denoted by light blue 
oval. Decreasing the tunneling amplitude between 72 and 73 and si- 
multaneously increasing the tunneling amplitude between 71 and 72 
can effectively transfer MF 71 — > 73. 

triangular loop that we will define as the product of tunnelings 
around the loop. Finally, we would like to note, that while the 
motivation of exchanging MFs is to be able to manipulate the 
information contained in topological qubits constructed from 
MFs in an effort to perform TQC, it is well-known that braid- 
ing by itself is insufficient for TQC. 6 However, MF exchanges 
are still crucual as one of the most direct tests of non-Abelian 
statistics and probably also for any future TQC schemes using 
MFs. 



III. MF TRANSPORT 

To understand how MF transport in a system of nanowires 
can be induced by tunneling, consider the simple system of 
nanowires shown in Fig. [Ha) consisting of three semicon- 
ducting nanowire segments. Two of these segments (shown 
by solid blue lines in the schematic in Fig. Q2b)) are in the 
TS phase, while the wire shown with the red dashed line is 
in the NTS phase and serves as the tunnel barrier connect- 
ing MFs. Each end of the wires in the TS phase supports a 
MF (shown as discs). In the initial state (shown in the upper 
panel of Fig. 02b)), the gate voltage of the NTS segment is 
chosen to allow a finite tunneling amplitude (shown by light 
blue oval) across it. This pairs up the MFs 72 and 73 into fi- 
nite energy states with a gap. Thus the operators 72,3 become 
gapped MFs (shown as dark blue discs) and cannot be used to 
store quantum information as can be done with the true zero- 
energy MFs (shown as light orange discs). The transfer of the 
MF from position 1 to 3 is achieved by adiabatically deacti- 
vating the tunneling in the NTS segment 2 — 3 and activating 
the tunneling in the TS segment 1 — 2. 

The process of transferring the MF from position 1 to 3 



shown in Fig. |TJb) is described by the time-dependent tun- 
neling Hamiltonian that is derived by extending the tunneling 
Hamiltonian in Eq.[T]and can be written as 

H = [Ci2a(t)7i72 + (23(1 ~ a(0)7273] (3) 

where (12 an ^ C23 are the activated tunneling amplitudes 
across the segments 1 — 2 and 2 — 3 respectively. Over the 
transfer process a(t) varies adiabatically from a(0) = 
to a(ti) = 1. It is convenient to understand the braid- 
ing procedure for MF operators in the Heisenberg repre- 
sentation 7j(i) = TP (t)jjU(t) where the U is the uni- 
tary time-evolution operator U(t) — Te~~ l $o H i T ) dT (which 
is a time-ordered exponential). The operators jj(t) can be 
computed from the Heisenberg equation of motion jj (t) — 
i[H^ HS> (i),7j(t)]. The Hamiltonian (Eq. |3]l describing the 
evolution of 7j(t) can be written compactly in terms of an 
effective B-field (B 3 (t)) as 

H {H) {t)= Yl tabcB a {t) lb {t) lc (t) (4) 
a,b,c=l,2,3 

where e a b c is the anti-symmetric Levi-Civita tensor. The time- 
dependent S-field given by 

B(t) = (l-a(i))C2 3 (l,0,0) + a(i)Ci,2(0,0,l). (5) 

The Heisenberg equation of motion for -f a (t) takes the form 

7a = 2e abc B b (t) lc (t). (6) 

This equation of motion is identical to that of the spin opera- 
tors a a (t) of a spin- 1/2 particle in a time-dependent magnetic 
field B(t) (with a Hamiltonian H^ H \t) = -B(t) ■ alt)). 
Furthermore, the initial condition on the operator ^(t) = 71 
corresponds to the spin-operator a(t) — cri(O) in an initial 
effective magnetic field B(0) = £23(1,0,0) that is aligned 
or anti-aligned with o\. Thus, after a time-evolution under 
an adiabatically varying magnetic field, the spin (and corre- 
spondingly the MF) remains aligned or anti-aligned with the 
final magnetic field B(t\) — £12(0, 0, 1) at time t — t\. This 
leads to the expression 

73(^1) = sgn(Ci 2 C23)7i(0)- (7) 

Thus the transfer of the tunneling amplitude from the segment 
2 — 3 to the segment 1 — 2, leads to transport of the MF from 
position 1 to position 3. The hopping of MFs between sites de- 
scribed by Eq.|7]is identical to the motion of regular fermionic 
operators under that action of tunneling. We will represent this 
process by the MF trajectory 

1 A 3. (8) 

The result in Eq. [7] is consistent with a somewhat different 
approach suggested by Kitaev. — 

IV. MFS AS DEFECTS IN DIMER LATTICES 

In the previous section, we saw that to transport a single 
MF from one position to another it was necessary to make use 
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FIG. 2: (Color online)(a) Configuration of TS nanowires in a square 
dimer lattice with 2 isolated MFs. MFs are effectively bound to de- 
fects (unpaired sites) on the dimer lattice. Tunneling between differ- 
ent TS segments fuses TS wires into isolated effective 'topological 
wires' shown in green (gray) blocks with MF end modes, (b) MF 
transfer processes analogous to Fig. Q] can extend effective TS wires 
and move MFs in dimer lattice in a way analogous to Alicea et al.— . 
(c) and (d): Similar transfer processes can be used to switch the MF 
end modes between different effective 'topological wires'. All op- 
erations of the MF square lattice either exchange, contract or switch 
ends of different 'topological wires'. 



of another pair of MFs which were coupled by a weak tunnel- 
ing so that they were not strictly zero-energy MFs. The single 
MF being exchanged had to be a true zero-energy MF with 
no tunneling, while the pair of MFs with tunneling between 
them may be thought of as a gapped MF dimer. Exchanges 
of true zero-energy MFs requires a generalization of this pic- 
ture to include several isolated MFs which are not coupled 
to any other MF by tunneling. As is clear from Fig. 02b), 
such a process also requires a supply of gapped pairs of MFs 
(i.e. MF dimers). In this paper, we will consider a system of 
nanowires with end MFs most of which are paired up by tun- 
neling into MF dimers as shown in Fig.[2fa). If all MF sites are 
completely paired up, then the system has no true zero-energy 
MFs and no topological degeneracy or non-Abelian statistics. 
Therefore we consider a system, where in addition to the MF 
dimers, there are a few isolated MFs that are unpaired by tun- 
neling. If one considers a regular lattice of nanowires so the 
the end MFs live on the vertices of a square lattice (Fig. [2]), a 
system of dimerized MFs forms a dimer covering of the lat- 
tice, while isolated MFs are associated with defects (unpaired 
sites) in the MF dimer lattice. As seen by comparing Fig. [2ja) 
and (b), the transport of MFs on the dimer-lattice is obtained 
by changing the dimerization pattern on the dimer lattice anal- 
ogous to Fig.QJb). 

The MFs on the dimer lattice can also be thought of as the 
end points of topological wires. To see this we define effective 
'topological wires', shown by green (gray) boxes in Fig. [2] by 
considering TS wire segments coupled by tunneling as a sin- 
gle topological wire. This identification, relates our proposal 
in a direct way to the proposal of Alicea et al. 26 However, the 
details of the physical implementation remain different and 
the dimer implementation presented in this paper will allow us 



to directly use Eq.|7] to determine the form of the braid-matrix 
in Eq. [2] The continuous processes required by Alicea et. al. 
for exchanging MFs were extending and contracting topolog- 
ical wire segments together with an operation that we will re- 
fer to as exchanging the ends of different topological wires. 
This process required bringing together a pair of topological 
wires in a tri-junction and effectively takes a pair of topologi- 
cal wires with end points 712 and 73 4 and creates a new pair 
of wires with end points 713 and 72.4. All these processes 
can be accomplished in MF dimer lattice by repeated applica- 
tion of the process shown in Fig.[TJassociated with Eq.|7] The 
analogue of the extension and contraction process in a dimer 
lattice is shown in Fig. |2ja) and (b), while the end switching 
process is shown in Fig. EJc) and (d). 



V. NON-ABELIAN STATISTICS OF MFS IN NANOWIRES 

In this section, we show explicitly, that exchange of MFs 
in any dimer lattice can always be described by an equation 
of the form of Eq. [2] Unpaired MFs can be exchanged via 
discrete tunneling operations of the form shown in Fig.[TJb). 
Since the physical positions of the MFs are exchanged by the 
correct sequence of MF transfers, the resulting transformation 
of the MFs at the end of the transformation t = t fi na i has the 
general form 



7l (t final) 
72 (t final) = 



= A 72 (0) 
A 7 i(0). 



(9) 



However, consistency with non-Abelian statistics also require 
us to prove that A A = —1. If A A = —1, the exchange trans- 
formation can be represented by the operator U of the form 
Eq.|H 

To show A A = —1, let us label the unpaired MFs, which are 
to be exchanged as 1 and 2 and take all other MFs as paired, 
(2n — 1, 2n) for n = 2,...,N, The positions of the MFs 
following each step (labelled by the index p) of the exchange 
process, which permutes the positions of the MFs, can be rep- 
resented by the function n p (j), where j = 1, . . . , 2N is the 
MF index. After each step p, the MF coordinates are updated 
from 7T p _ 1 to TT p according to the relation 



(10) 



where C p is a cyclic (clockwise or anticlockwise) permutation 
of the MFs a p , 2n p — 1 and 2n p corresponding to Eq.|7] Here 
we choose a p to be one of the unpaired MFs, 1 or 2, and n p > 
1 such that MF dimer (2n p — 1, 2n p ) is paired. The equation 
of motion for the unpaired Majorana operators corresponding 
to Eq.|7]is 



J7t p+1 (a p ){t p +l) - \^ P (a p ){t P )- 



(ID 



where 



Ap — Sgn(C 7rp+1 (2n p -l) 7 r p+1 (2n p )CTr p (2np-l)Tr p (2n p ))- (12) 

The total sign AA picked up by the unpaired MFs is the prod- 
uct A A = Yl p Ap- To calculate this product we define a se- 
quence Q p = sgn(n„>i C7r p (2n-i), 7 r p (2n))- From Eq. [TJ it 



follows that Q p +i 



-- XpQp so that 

Q final = X\Q p= Q. 



(13) 



Note that since each cyclic permutation C p contains even 
number of exchanges (i.e is an even permutation), the per- 
mutations 7T p at each step (including the final permutation 
Kfinai), which is a product of C p s, is also an even permu- 
tation. 

Since the Hamiltonian is required to return to its initial 
configuration, the MFs at positions (2n — 1, 2n) must be 
paired by tunneling for n > 1. This requires that tt final is 
composed of a pair exchange of the positions of MF dimers 
(2n — 1, 2n) o (2n' — 1, 2n') together with possible internal 
flips (2n — 1 •<-> 2n) of the dimers. Since, n final is an even 
permutation, and dimer exchanges are even permutations, the 
number of internally flipped dimers (2n - 1 f> 2n) in it final 
is even. Moreover, the unpaired pair of MFs (1,2) is flipped 
in it final- Thus an odd number of the paired MF dimers 
(2n — 1 ■<-> 2n), must be flipped for n > 1. Each such dimer 
flip changes the sign of Qf 

inali since C,2n— l,2n — C2n,2n— 1 

for n > 1. This leads to the relation Q final = —Q P =o — 
XXQ p= q, proving the consistency condition for non-Abelian 
statistics i.e. AA = — 1. 



VI. EXCHANGE AROUND A TRIANGULAR LOOP 

A specific realization of the non-Abelian statistics in 
nanowire systems is provided by a triangular loop geometry 
shown in Figs. [3] and [4] The triangular loop consists of one 
end (A 2 ,B 2 and C 2 ) of each of three TS segments (A, B and 
C) connected by NTS segments to form a triangle. The other 
ends are labeled A\, B\ and C\. The MFs to be exchanged, 
referred as 1 and 2, are assumed to be localized at 2 of these 6 
ends of TS segments. Each of the steps for the MF exchange 
(shown in Figs.[3]and |4]) consists of moving exactly one MF 
from one position to the other (shown by dotted arrows) by 
adiabatically turning off the tunneling in some wire segment 
and increasing it in an adjoining segment as discussed before. 

The procedure to exchange the MFs 1 and 2 at the ends 
of different TS segment through the tri-junction takes place 
in four steps shown in Fig. [3] The signs associated with the 
exchange A and A can be determined by following the trajec- 
tories of the MFs 1 and 2 and applying Eqs.Qand [9] From 
Fig. |3] it is clear that the sequence of positions followed by 
the MFs 1 and 2 are 



MF 1: A 2 ■% B 2 

(3) 
C 2l 



MF 2: B 2 ^Ci A 2 

(2) (4) = (1) 



(14) 



respectively. Here we show only the MF that is moved in each 
step, which is numbered in Fig. |4]as (j — 1, . . . , 4) (marked 
below the arrows in EqlPfli. The MF motion is shown using 
the notation defined in Eq.|8]so that the sign can be calculated 
using Eq. [7] Applying Eq. [7] the parameters A and A simplify 
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FIG. 3: (Color online) MFs 1 and 2 at the ends of different TS seg- 
ments are exchanged. This is achieved by switching tunnelings on 
and off on TS and NTS segments in 4 steps going from a state shown 
in one panel to the next panel. Dotted arrow shows motion of MF 
from the previous panel. The labelling for the sites Ai,A2, B\,B-i, 
G\ and C2 is shown in panel (1). 
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FIG. 4: (Color online) MFs 1 and 2 at the ends of the TS segment 
on the left leg are exchanged in seven steps similar to Fig. [3] Step 
(7) transfers state shown in panel (6) back to panel (1) with the effect 
that the Majoranas 1 and 2 are interchanged. 



to 



A = -A = sgn(C A2 s 2 )x 



(15) 



where x = sgn(CA 2 s 2 C-B2C 2 Cc 2 A 2 ) is defined to be the chi- 
rality of the tri-junctioniS 

Similarly MFs at the ends of the same TS segment can be 
exchanged using six steps shown in Fig. [4] From Fig. [4] it is 
clear that the sequence of positions followed by the MFs 1 and 
2 are 



MF 1: Ai 
MF 2: A 2 



A 2} 


c 2 


B 2 




B 2 


(3) 




(4) 




(7)^(1) 


C 2> 




C 2> 


B 2 


^A 


(2) 




(5) 




(6) 



Ao 



(16) 



respectively. The step (7) is not explicitly shown in Fig. |U 
since it is equivalent to (1). Applying Eq.|9] the parameters A 
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and A simplify to 

A = -A = sgn[( Al A 2 )x (17) 

where \ is the junction chirality. 

Thus, using Eq.[T7]and Eq.|9] we obtain the the result that 
the unitary time-evolution of the MFs 71 and 72 under ex- 
change can be described by the unique braid-matrix 

JJ = gf Xsgn(Cl2)7l72 (Jg) 

where C12 is the tunneling amplitude of the segment separating 
71 and 72. The quantities (12 and x for a specific network are 
calculated in the appendix. 



VII. CONCLUSION 

Non-Abelian statistics for MFs at the ends of TS nanowire 
segments can be realized by introducing time-varying gate- 
controllable tunnelings between MFs in a nanowire system to 
exchange the end MFs. Similar to the previous proposal for 
braiding MFs in ID wires, our system can also be embedded 
in 3D leading to the possibility of non-Abelian statistics in 
3D. The isolated MFs being exchanged in the tunneling ge- 
ometry considered in this paper may be thought of as defects 
in a dimer lattice i.e. sites that are unpaired by tunneling. 
Alternatively, this system may also be thought of as a dis- 
cretized implementation of the continuous nanowire network 
proposal of Alicea et. ali2& However, the discrete implemen- 
tation discussed in this paper allows us to compute the braid 
matrix explicitly in terms of MF overlaps. The non-Abelian 
braid matrix for exchange around a triangular loop geometry 
is given by a product of the fusion channel of the MFs (1.2 
and the junction-chirality x- The fusion channel (12 is sim- 
ply the tunneling matrix-element between the MFs being ex- 
changed and the junction-chirality is the product of tunneling 
terms around the triangular junctions. Thus the braid-matrix 
in the tunneling geometry considered in this paper is com- 
pletely determined in terms of microscopic tunneling param- 
eters by Eq. [18] making nanowire systems a well-controlled 
platform to realize non-Abelian statistics. 
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FIG. 5: (Color online)Schematic of orthogonal nanowire system on 
a superconductor (shown as light rectangle) that generates tunneling 
of MFs (shown as light orange discs). The entire system is subject to 
an in-plane magnetic field to generate Zeeman coupling. Nanowire 
segments in the TS phase are shown as dark blue rectangles with 
end MFs. Tunneling is generated between MF T and MF L by con- 
ventional tunneling across a nearly depleted nanowire in the NTS 
phase. The tunneling can be calculated using the Bardeen tunneling 
formula 3 - as the matrix element of the current operator in the mid- 
dle of the wire (black dotted line). Similarly tunneling is generated 
between MF L and MF R by lowering the topological gap so that 
the wave-functions have significant overlap at the middle of the wire 
(black dotted line). 

Appendix A: Calculation of tunneling matrix elements for a 
specific nanowire system 

In the main text of the paper we saw how controlling the 
tunneling between MFs can be used to generate transport of 
Majorana fermions from one point to another and eventually 
generate exchanges and braids that are useful for TQC. The 
sign of the resulting exchange was found to be determined by 
the signs of various tunneling matrix elements. While the ex- 
istence of non-Abelian statistics is demonstrated in the paper 
in general, the signs of the tunneling matrix elements them- 
selves are depend on the microscopic details of the system. 

In this appendix, we calculate the tunneling matrix ele- 
ments between the various MFs at a junction for a network 
of orthogonal wires on a superconducting substrate as shown 
in Fig. [5] A Zeeman potential is applied at 45 degrees to the 
wires. The wires in the TS phase (shown in dark blue), which 
support MFs (shown as orange discs) at their ends, are taken 
to have a Rashba spin-orbit coupling generated from interac- 
tion with the superconducting substrate. The Bogoliubov de 
Gennes (BdG) Hamiltonian for the wire along x is given by 

H B dG = {~vdl~Kx))T z +V z cr-B+iad x a y T z +AT x (Al) 
and for the wire along y is given by 

H B dG = {-vdy - v(y))r z + V z a- ■ B - iad y cr x T z + At x . 

(A2) 

The direction of the Zeeman field is B = (x+y) /\/2. Follow- 
ing the spin-rotation and phase transformations in Ref. 1I20I1 . 
for negative Rashba coupling a < 0, the Majorana wave- 
functions at the left and the right ends of the x wires have 
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the form 



ia y u{x)e- l ' t '' 2 



and d>p = ia„ 



KJyui-xy^/ 2 
(A3) 



respectively, where (b — sin -1 and u(x) is a real 2- 
spinor. Note that in this geometry there is now an additional 
condition for the wire to be gapped i.e -\/A 2 + /i 2 < Vz < 
A\/2. This constraint implies j < (b < -|. The Majorana 
wave-functions for the Majorana fermions at the bottom and 
top ends of the wires parallel to the y axis have the form 



Q 



-icp/2 



u{y)e 
i(j y u{y)e^ 2 



and 



iQa x 



u{y)e 
ia y u(y)e~ 



(A4) 



respectively, where Q = e ~ l7rCT -/ 4 . 

Transport of MFs is generated by introducing tunneling into 
the system of MFs shown in Fig. [5] The junction chirality 
X defined in the paper depends only on the tunneling from 
3 end MFs MF {L, R, T} with wave-functions 4>l,r,t- Let 
us start by considering the MF overlap across the NTS seg- 
ments (shown as red dashed lines in Fig. [5} which is simplest 
to understand in the limit of low negative chemical potential 
fi = where Vz, A. Physically, this corresponds 

to a wire that is nearly depleted of electrons and only acts as 
a tunnel barrier. In such a case, the Majorana wave-function 
in the NTS wire has the usual exponentially decaying form 
^(x) = ^(a^e -7 ' 1-2 ^' as in a barrier, where 7 ~ y / 2m\fj,\ 
and xj is the position of the interface between the TS and 
NTS wire segments. The tunneling matrix elements Qj be- 
tween two MFs at x/.i = —a/2 and xi >2 — a/2 across the 
NTS wire can be calculated from the matrix-elements of the 
current operator and the wave-functions in the middle of the 
wire^ (x = as shown by the dark dotted lines in Fig.|5]l 

C = ^\(0)T z d x y 2 (x)\ x =o - d x y\(x)\ x = T z y 2 (0)} 



-^(-2)^1(2) 



(A5) 



where p 



-7a 



is the overall tunneling strength and we 



have assumed A 3> a. In this limit, the overlap between a pair 
of Majorana wave-functions $1 = (ui(x), ia y u\ (x)) T and 

^2 = (u2{x) , ta y U2(x)) T is given by M = 2iplm(u\u2)- 
This is purely imaginary and manifestly anti-symmetric as 
expected. Furthermore since the fundamental spinor u(x) in 
terms of which each of u\ 2 are written is real, we can write 



it as u = (cos#,sin(9) T , where the parameter 6 depends on 
Vz, /i, a etc. With the help of these relations it is easy to tab- 
ulate the Majorana tunneling matrix as an anti-symmetric ma- 
trix for the states in the order (L, R, T) as 



'V2 



y/2 cos <p sin 29 sin 20 

* (sin (b + cos (b cos 26) 

* * 



(A6) 

where the elements in the * have been left empty since they 
are determined by the anti-symmetry constraint. The junc- 
tion chirality \ m me previous section, used only the Majo- 
rana modes L,T,R and is calculated using the expression, 
X = CrlCltCtr = p 3 cos 2 (b sin 2 2#[cos 28 + tan </>], which 
is always positive, since tan</> > 1 for the Zeeman direc- 
tion B = (x + y)/V2 and negative Rashba coupling a < 0. 
Changing the Zeeman potential to B = (x — y)/v2 flips 
the chirality. Changing the sign of the Rashba coupling a 
requires us to change <pL — > i(u(x)e l ^ > ^ 2 , — io y u{x)e~ 1 '^ 2 ). 
Since all the other wave-functions are derived from symmetry 
transformations applied to (p^, the rest of the calculation goes 
through as is with the only difference that u{x) changes to 
iu(x). Therefore the final result for the chirality of the junc- 
tion is independent of the Rashba coupling. 

The signs acquired by MFs on exchange is dependent 
also on the tunneling between the MFs MF L,R across a 
TS segment. The tunneling amplitude between Majorana 
fermions on the same topological segment is well-controlled 
and can be calculated in the limit of a long topological 
wire (wire length L > a/Vz)- In this limit, the MFs 
only overlap in the limit where the gate potential is tuned 
so that the wire is driven towards a phase transition by 
tuning p near to /i = \fV% — A 2 . The relevant slow- 
est decaying spinor component then determines the tun- 
neling matrix elements and is given by u(x) — {Vz + 



S gn{a)JVj^-ti T exp{-^{JvJ=p? - A)). 



The 



tunneling matrix element is given by M <~ — i(x>S> ^ayT^ r = 
— iaRe[u L <r y Ufi\. Substituting u, we find the over- 
lap to simplify to M oc ~iacos(f)e~ 2L '€ where £ = 
|a|/(-^/V| — A 2 — /i). Thus the sign of C,l.r is determined 
by the sign of the Rashba spin-orbit coupling a. Using these 
results together with Eq. lA6l one can check the physically rea- 
sonable result that changing the sign of the Rashba coupling 
a, flips the sign of the MFs on interchange. 

Thus, provided care is taken to ensure the conditions for 
the braid discussed in this section, the sign of both clock- wise 
exchanges is positive for positive Rashba coupling a and neg- 
ative other-wise. For a given Rasbha coupling the sign of the 
braid can be altered by considering anti-clockwise braids. 
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